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^ ' Abstract. The kinematic structure of a sample of planetary nebulae, consisting of 23 [WR] central stars, 21 weak emission 

ON , line stars (wels) and 57 non-emission line central stars, is studied. The [WR] stars are shown to be surrounded by turbulent 

nebulae, a characteristic shared by some wels but almost completely absent from the non-emission line stars. The fraction of 
objects showing turbulence for non-emission-line stars, wels and [WR] stars is 7%, 24% and 91%, respectively. The [WR] stars 
show a distinct IRAS 12-micron excess, indicative of small dust grains, which is not found for wels. The [WR]-star nebulae are 
on average more centrally condensed than those of other stars. On the age-temperature diagram, the wels are located on tracks 
of both high and low stellar mass, while [WR] stars trace a narrow range of intermediate masses. Emission-line stars are not 
found on the cooling track. One group of wels may form a sequence wels— [WO] stars with increasing temperature. For the other 
' groups both the wels and the [WR] stars appear to represent several, independent evolutionary tracks. We find a discontinuity 

in the [WR] stellar temperature distribution and suggest different evolutionary sequences above and below the temperature 
gap. One group of cool [WR] stars has no counterpart among any other group of PNe and may represent binary evolution. A 
prime factor distinguishing wels and [WR] stars appears to be stellar luminosity. We find no evidence for an increase of nebular 
"qJj expansion velocity with time. 
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1. Introduction ary sequence [WC1 1]->[W02] (Zijlstra et al. 1994, Acker 

, „ , , , _ T s , et al. 1996, Pena et al. 2001), followed by PG1159 stars 

Some central stars (cores) or planetary nebulae (PNe) show .... . , . , „ XT . 

. . ,. ... „ 7 ,„„ (Werner et al. 1992). Acker & Neiner (2003) propose a se- 

broad, stellar emission lines similar to Wolr-Rayet stars. mfrl1l1 ,,,,„., mi^Ai mmn i-> ^ 

• . . J . quence: [WC1 1]^[WC4H[W04]->[W01]. Parthasarathy et 

Although superficially similar, they diner from classical WR . ,,„„„, , . , . , , , , , 

° r J J al. (1998) have suggested that the wels are related to the 

stars in their degenerate structure, much lower masses, a wider „„,,.„ . ru ,™ , r ,,,„ rl , m „. . „„,,„ „ , 

b . . i-i PG1 159 stars: [WCL]-^[WCE]-*wM/PGl 159—>PG1 159. But 



range of temperatures, and being limited almost exclusively to 
carbon-rich stars (there are very few counterparts to the mas- 
sive WN stars). Similar to their massive WC-star counterparts, 
they are deficient in hydrogen. They form a class named [WC]- 
type which historically was further divided into early [WCE] 
and late [WCL] groups. A less numerous class of PNe central 
stars show narrower and weaker emission lines than the [WR]- 
type stars. These are named wels (weak emission-line stars), as 
defined in Tylenda et al. (1993). 

The evolutionary status of the [WR]-type stars is still very 
uncertain, and it is unclear whether there is any evolution- 
ary relation to the wels. Nebular data suggest an evolution- 



this sequence is by no means proved; the precise location 
of the wels in relation to the [WR] stars is under discussion 
(Marcolino & de Araujo 2003). Pena et al. (2001) also ar- 
gue against too closely identifying wels with PG1159 stars 
(see also Koesterke 2001). Not all PG1159 are hydrogen-poor 
(Dreizler et al. 1996), showing that not all PG1159 stars will 
have evolved from [WR] stars. Neither are all wels hydrogen- 
poor: a dual evolutionary sequence may be expected. 

The lack of hydrogen is often taken to indicate that the 
[WR] stars have undergone a late thermal pulse (or helium 
flash), either during the post-AGB evolution (a so-called Late 
Thermal Pulse or LTP) or on the white dwarf cooling track 
Send offprint requests to: K. Gesicki (Very Late Thermal Pulse or VLTP: Herwig 2001, Hajduk et 
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al. 2005). Following such a pulse, the star rejuvenates and re- 
traces part of its earlier evolution. This predicts that nebulae 
around [WR] stars should be more evolved than around non- 
[WR] stars, but this is not confirmed by observations; the prop- 
erties of the PNe around [WR] stars do not differ from those 
around other central stars (Gorny 2001). 

In this paper we discuss a much larger sample of [WR] stars 
and wels than has previously been available. We apply velocity- 
field analysis to locate the objects on derivatives of the HR di- 
agram. In Sect. 2 we describe the methods applied for analysis 
of PNe with specific attention towards automatic model fitting. 
Sect. 3 presents the 33 newly analyzed PNe. In Sect. 4 we dis- 
cuss the nebular and stellar parameters for the full sample of 
101 PNe. In Sect. 5 we discuss possible evolutionary relations. 

Crowther et al. (1998) have refined the WC and the WO 
schemes used for WR stars and define a unified classification 
for massive WR and low-mass [WR] stars. Acker & Neiner 
(2003) developed this scheme for a sample of 42 PN central 
stars, classified into [WO 1-4] and [WC 4-11] stars. This clas- 
sification (used in the present paper) is based on the ionization 
level of the elements (depending on the wind-temperature of 
the PN core), showing essentially carbon lines for the coolest 
stars and oxygen lines for the hottest ones. 

2. Methods 

We derive information on the nebulae using a combination of 
line ratios, diameters, and high resolution spectra. The diam- 
eters and line ratios are used to fit a photo-ionization model. 
The model is constrained to reproduce the intensity distribu- 
tion of images, if any are available. Line profiles are obtained 
from the high resolution spectra, and are fitted using the emis- 
sivity distributions of the photoionization model, and assuming 
a velocity field. This is done using the Torun models (Gesicki 
et al. 1996, 2003). The velocity fields include separate contri- 
butions from expansion and turbulence. Turbulence always in- 
cludes the instrumental broadening and the thermal broadening 
(calculated from the photoionization model), but some objects 
may require additional turbulence. 

The models assume spherical symmetry. Some strongly 
bipolar objects cannot be fitted, because they show very irreg- 
ular velocity profiles. But the majority of objects can be well 
reproduced with a spherical model. The model corrects for the 
size of the aperture used for the spectroscopy, and includes see- 
ing effects. 

The diameters and distances are adopted from the litera- 
ture. The models find the density distribution, stellar temper- 
ature and the velocity field. The density as function of radius 
is usually assumed in the shape of a reverted parabola, how- 
ever if images were available we compared the computed sur- 
face brightness and improved the run of the density. The stellar 
temperature assumes a black-body spectrum energy distribu- 
tion - therefore we prefer to call it 7Vb instead of r e ff. This 
assumption can be challenged, especially for the [WR] cen- 
tral stars where the comparison between the stellar temperature 
T, obtained from non-LTE modelling and T2/3 which refers to 
the radius R(t Ross = 2/3) shows differences. The differences 
are significant for hotter [WO] stars (Koesterke & Hamann 



1997) while smaller for cooler [WC] objects (Leuenhagen et 
al. 1996). The general tendency is always the same: the T2/3, 
T e g or Zanstra temperatures are smaller than T„. Despite this 
discrepancy we apply the 7Vb for further discussion since it 
has the advantage of providing a uniform measure of the tem- 
perature over a range of stellar classes. 

The velocity distribution is assumed to vary arbitrarily and 
smoothly with radius. In Gesicki & Zijlstra (2000) we com- 
pared the Torun model analysis with more traditional methods 
of deriving the expansion velocities from the spectra. Since 
long time it was known that different observed lines resulted 
in different velocities (e.g. Weinberger 1989) indicating a ve- 
locity gradient. Our model has this advantage that it combines 
these data into a single velocity curve. Earlier it was not ob- 
vious but the Torun models have shown that even the shape 
of a single line can indicate a velocity gradient. Good veloc- 
ity fields can be obtained if a larger number of lines are avail- 
able. The resulting fields show detailed structure, with velocity 
peaks at the outer and sometimes at the inner radius of the neb- 
ula (Gesicki & Zijlstra 2003). If fewer lines are available, the 
velocity field is less constrained and simplifying assumptions 
need to be made. From the full velocity curve V(r) we derive a 
single parameter characterizing the nebular expansion. We de- 
fine the expansion velocity as a mass-weighted average over 
the nebula, V av . This parameter has been shown to be robust 
against the simplifying assumptions: it can be accurately deter- 
mined even when the velocity field itself is uncertain (Gesicki 
et al 2003). This allows us to define a kinematic age to the neb- 
ula. 

2.1. The genetic algorithm 

The Torun models were recently adopted to automate the 
search for a best-fit model, using an optimization routine. We 
applied the method widely known as the genetic algorithm. 
While it is still not a very popular optimization technique, it 
has proved to be very effective and robust in many problems, 
in particular of astrophysical origin. For a review, we refer to 
a paper by Charbonneau (1995) who is the author of the pub- 
licly available code PIKAIA for a genetic algorithm 1 , used for 
computations in this work. 

The genetic algorithms have been invented as an optimiza- 
tion technique mimicking the processes of biological evolution 
(e.g. Koza 1992). They lead to the selection and adaptation 
of life forms to the conditions of the natural environment. In 
this sense the evolution can be thought of as a powerful opti- 
mization algorithm. The genetic optimization starts with a set 
(population) of randomly chosen individuals (parameters of the 
mathematical model encoded in the 'genomes'). The whole pa- 
rameter space of the problem creates the environment for these 
individuals. The likelihood of survival of a given individual 
is determined by its 'fitness' function /; usually, in the least 
squares minimization, / = 1/ >Jx*. Then the fitness of a par- 
ticular individual is a measure of goodness of fit to the studied 
data set. The population evolves through a number of gener- 
ations. At every generation, the genetic algorithm evaluates / 

1 http://www.hao.ucar.edu/Public/models/pikaia/pikaia.html 
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Fig. 1. The best fit luminosity and temperature for photoioniza- 
tion models for a single object (He2-113, PNG 321.0+03.9) 
for 40 different PIKAIA runs. The larger the size of the circle 
the better the fit quality. 

resulting from each parameter set and processes information 
by applying to the genomes genetic operators like crossover, 
mutation and selection. The best fitted members of the pop- 
ulation are used to produce a new generation and the process 
continues until some convergence criteria is reached. Thanks to 
keeping the memory of the best fitted individuals, the genetic 
algorithms are much more efficient than Monte Carlo based 
search techniques. They are robust and especially well suited 
for multi-dimensional problems, models possessing multiple 
local extrema and discontinuities. They are non-gradient meth- 
ods and thus especially well suited for models dependent in a 
sophisticated way on their parameters. Details of practical ap- 
plication of the technique are given in an excellent introduction 
by Charbonneau (2002). 

The genetic optimization replaced the 'trial and error' 
search of parameters applied in our earlier publications (e.g. 
Gesicki & Zijlstra 2003, Gesicki et al. 2003). The algorithm 
appeared very effective in comparison with the earlier work. 
We note also that because of random initialization of the pop- 
ulation the method is now free from initial biases. The search 
procedure again follows two steps. First PIKAIA searches for 
the values of the nebular mass and the stellar temperature and 
luminosity which optimize the fit to the observables, before at- 
tempting to find the velocity field which optimizes the fit to the 
line profiles. The following expression is optimized: 
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where 0(f) are N observations (i = 1 , . . . , N), C(i) are model 
predictions dependent on p parameters and cr(i) are individ- 
ual errors of the observations. In some cases considered in this 



work, N is very small and the determination of uncertainties on 
the fitted is difficult (and without much meaning in terms of the 
formal, statistical approach). Therefore we did not perform a 
formal error analysis of the results. Nevertheless, in a well per- 
formed optimization search the error analysis is as important as 
the determination of the fit parameters. 

To check the possible degeneracy of the solutions and to ob- 
tain an idea about the errors, we collected the values to which 
the search process converged in many independent runs with 
the same start, and projected them on chosen parameter planes. 
The distribution of these parameters gives insight on the sig- 
nificance of the minima, whether they can be well localized 
and fixed in the parameter space of the problem. In Fig.^we 
present an example of fitting the photoionization model for a 
single PN (He 2-1 13, PN G 321.0+03.9). Results are plotted for 
40 different runs of PIKAIA, using the same set of observables. 
A scatter in the obtained parameters is seen: few of the points 
are overlapping. The size of the symbols corresponds to the fit 
quality. For this PN the temperature is better constrained than 
the stellar luminosity. The situation for other PNe is similar. For 
higher temperatures the models seem to be better constrained. 
We estimated the accuracy in 7Vb as +2000K and in log L/L 
as +0.3. 

Similar test were performed for the velocity procedure. 
When fitting a simple velocity field (linear or parabola-like) 
the PIKAIA routine repeats very well the results. More com- 
plicated velocity fields result in different shapes for different 
runs; however the largest discrepancies appear in the least con- 
strained parts of the velocity curves. Nevertheless the mass- 
averaged expansion velocity remains well determined and is 
similar in all runs. From the distribution of models, we esti- 
mate the accuracy of V av as about ±2km s . 



2.2. Observations 

A large number of PNe have been analyzed in previous pa- 
pers using the Torun models. In most cases the velocity fields 
were determined from three lines (hydrogen Ha 6563 A, [N n] 
6583A and [O m] 5007 A) or fewer. To this sample we add in the 
present paper 27 unpublished CAT observations, which cover 
only Ha and [N n] 6583A, supplemented with recent echelle 
NTT observations of 6 PNe with emission-line central stars. 
PNe with new data are listed in Table 1 . 

The ESO Coude Auxiliary Telescope (CAT) was a 
subsidiary 1.4m telescope feeding into the Coude Echelle 
Spectrograph (CES) located at the neighboring 3.6m telescope. 
The CAT observations of 27 objects were performed during 
1993 and 1994. The long camera was used, giving a spectral 
resolution of 60 000 (corresponding to 5 km s _1 ). The slit width 
was 2". The observations use a long slit, sampled with 2" pix- 
els. However, the spectra analyzed here use only the central 
row of pixels. (As the nebulae studied here are compact and 
the CAT was not an imaging-quality telescope, no spatial in- 
formation was expected.) The spectrum covers one order of the 
echelle, covering Ha and the [Nn] lines at 6548 A and 6853 A 
but no other nebular lines. 
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Fig. 2. The observed and modelled lines H I 6563 A and [N n] 
6583A. The circles correspond to the observed profile, the line 
to the fitted model. The line fluxes are normalized to unity. The 
X-axis velocity scale given in the lowest boxes is the same for 
all plots. 



Six objects with emission-line central stars were observed 
with the ESO New Technology Telescope (NTT) during June 
2001. The echelle spectra were obtained using ESO Multi 
Mode Instrument (EMMI), with grating 14 and cross disperser 
3, giving a resolution of 60000. The slit width was 1". The 
spectra covered the wavelength range 4300-8 100 A. The spec- 
tra were summed over the slit length of 3". Exposure times 
were typically 120 seconds. Between 3 and 7 lines per object 
showed high enough S/N to be used in the velocity analysis. 




-60 -30 30 60- 
[km/s] 



60 -30 30 60 
km/ s] 



Fig. 3. The observed and modelled lines Hi and [Nil]. 
Continuation of Fig. [2 



The details of the data reduction are given in Gesicki & Zijlstra 
(2003). 

The profile fitting procedure requires specification of the 
center of the emission line with zero velocity. We did not per- 
form the full radial velocity correction. Instead we assumed 
the zero position in the centre of symmetry for the almost sym- 
metric lines and usually located the zero position in the centre 
between the two maxima of the asymmetric lines. In the sec- 
ond case the presence of two (or more) spectral lines helped to 
locate the zero position, the obviously hopeless cases were not 
analyzed. 
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Table 1. The expansion velocities and other data concerning the nebulae and the central stars, based on observations described 
in this paper. The column 'remarks' indicate for PNe observed with EMMI, how many emission lines were used for the analysis. 



PNG Name log 7Vb logL Dist. R ollt Mj on V„ V tr b M co „ t t i y „ Central Remarks 







[K] 


[L Q ] 


[kpc] 


[pc] 


[Mo] 


[km s 


'] 


[M ] 


[kyrs] 


star 




000.9-04.8 


M 3-23 


5.18 


3.0 


4.0 


.11 


.10 


24 


13 


0.61 


5.1 






002.6-03.4 


M 1-37 


4.40 


3.9 


8.0 


.04 


.08 


27 





0.60 


1.7 


[WC11]? 




002.7-04.8 


M 1-42 


4.92 


3.0 


4.0 


.08 


.16 


13 





0.60 


5.4 






003.6+03.1 


M 2-14 


4.64 


3.0 


8.0 


.05 


.06 


17 


10 


0.60 


2.9 


wels 


6 lines 


004.6+06.0 


H 1-24 


4.56 


3.8 


7.0 


.08 


.13 


24 





0.58 


3.7 


wels 




006.4+02.0 


M 1-31 


4.76 


3.8 


8.0 


.06 


.29 


19 





0.61 


3.2 


wels 




008.2+06.8 


He 2-260 


4.80 


3.2 


12.0 


.06 


.13 


17 





0.61 


3.5 






211.2-03.5 


M 1-6 


4.60 


3.5 


4.0 


.03 


.05 


24 





0.62 


1.4 






217.4+02.0 


St 3-1 


4.96 


2.5 


4.0 


.14 


.20 


26 





0.60 


6.1 






232.0+05.7 


SaSt 2-3 


4.80 


2.5 


4.0 


.02 


.01 


19 





0.64 


1.1 






232.4-01.8 


M 1-13 


5.04 


2.9 


4.0 


.10 


.19 


14 





0.60 


6.5 






253.9+05.7 


M 3-6 


4.72 


4.0 


1.5 


.036 


.03 


18 





0.61 


2.0 


wels 




258.1-00.3 


He 2-9 


4.73 


2.9 


1.5 


.015 


.01 


25 





0.65 


0.7 


wels 




261.6+03.0 


He 2-15 


5.17 


3.0 


2.2 


.13 


.32 


22 





0.71 


6.4 






283.3+03.9 


He 2-50 


5.12 


2.9 


5.1 


.14 


.32 


17 





0.59 


8.1 






292.4+04.1 


PB 8 


4.76 


3.6 


5.0 


.06 


.18 


22 


9 


0.61 


2.9 


[WC5-6] 




300.7-02.0 


He 2-86 


4.83 


3.4 


3.0 


.03 


.05 


14 


1 1 


0.62 


2.0 


[WC4] 




304.8+05.1 


He 2-88 


4.72 


2.8 


4.0 


.03 


.03 


23 





0.62 


1.4 






309.0-04.2 


He 2-99 


4.49 


2.7 


4.0 


.16 


.25 


45 


10 


0.57 


4.6 


[WC9] 




316.1+08.4 


He 2-108 


4.51 


3.8 


4.0 


.10 


.14 


21 


9 


0.57 


5.0 


wels? 




321.0+03.9 


He 2-113 


4.48 


2.8 


2.0 


.01 


.003 


18 


15 


0.64 


0.6 


[WC10] 




325.0+03.2 


He 2-129 


4.87 


3.8 


7.0 


.03 


.08 


18 


13 


0.63 


1.7 






325.8+04.5 


He 2-128 


4.70 


3.3 


5.0 


.06 


.13 


12 


10 


0.59 


4.3 






347.4+05.8 


H 1-2 


4.96 


3.9 


7.0 


.02 


.06 


13 


10 


0.64 


1.4 


wels 


3 lines 


351.1+04.8 


M 1-19 


4.72 


3.2 


11.0 


.10 


.20 


26 





0.61 


4.3 


wels 


5 lines 


352.1+05.1 


M2-8 


5.11 


3.2 


7.0 


.06 


.16 


14 


14 


0.61 


3.9 


[W02-3] 


7 lines 


355.9+03.6 


H 1-9 


4.58 


4.0 


9.0 


.04 


.09 


20 





0.61 


2.1 






356.1+02.7 


Th 3-13 


5.05 


3.3 


8.0 


.03 


.05 


32 





0.65 


1.1 


wels 


3 lines 


356.5-02.3 


M 1-27 


4.41 


3.6 


3.0 


.04 


.05 


21 





0.60 


2.0 


[WC11]? 




357.1+03.6 


M3-7 


4.72 


3.2 


4.0 


.05 


.07 


23 





0.61 


2.4 


wels 


3 lines 


357.4-03.2 


M2-16 


4.95 


3.3 


7.0 


.09 


.25 


21 


14 


0.61 


4.5 






357.4-03.5 


M2-18 


4.64 


3.5 


7.0 


.06 


.09 


24 





0.60 


2.8 






359.1-02.3 


M3-16 


4.57 


3.8 


7.0 


.08 


.14 


30 





0.59 


3.1 







3. The modelled planetary nebulae 

Figures |2jf4]present the observed and modelled nebular emis- 
sion lines H I 6563A and [N n] 6583A. Although for the indi- 
cated in the last column of the Table^six PNe more lines were 
used for model fitting ([Om] 5007 A and in some cases [Oi] 
6302A, [S n] 6732A, [Sra] 63 11 A, Hen 4686A), we present 
only the two lines for a consistent presentation. We do not draw 
the model radial distributions of the velocity, density and sur- 
face brightness because our main results are the mass-averaged 
expansion velocities. The [WR]-subclass and wels classifica- 
tion is adopted from Acker & Neiner (2003) and for the low- 
est temperature [WC11] stars from Gorny et al. (2004). The 
summary of nebular and central star parameters are listed in 
Tabled The last column indicates whether more than two lines 
were used for modelling. A few objects are discussed in the 
following. 

PNG002.6-03.4 (M 1-37), [WC11]"?"). The density 
structure is adopted to agree approximately with the image of 
Sahai (2000). The brightest nebular structure is the inner ellip- 



soidal ring 1.5 x 2.5" which is modelled by our 2" sphere. The 
multipolar lobes as well as the extended spherical halo are be- 
yond our analysis. We obtained a lower temperature than Zhang 
& Kwok (1991) for this low excitation nebula. 

PNG253.9+05.7 (M 3-6, wels). The radio image (Zijlstra 
et al. 1989) reveals a complicated structure with four maxima. 
Our spherical model is density bounded, but the line ratios are 
not well fitted and suggest a mixture of ionization and density 
boundaries. A small central cavity is indicated. The best ve- 
locity field which reproduces the multicomponent structure of 
the [N n] line shows multiple maxima. Considering the compli- 
cated nebular image our model should be treated with caution, 
but because the spectral lines are almost symmetric the mass- 
averaged expansion velocity is believed to be reliable. 

PNG321.0+03.9 (He 2-113, [WC 10]). The HST image 
(Sahai et al. 2000) reveals complex, highly aspherical struc- 
tures which makes spherical modelling doubtful. Nevertheless 
the average velocity still is a useful parameter. In our observa- 
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Fig. 4. The observed and modelled lines Hi and [Nil]. 
Continuation of Fig. [3] 

tions this PN is unresolved and the emission lines are symmet- 
ric. 

4. Derived parameters 

The newly analyzed data is combined with an earlier sample 
discussed in Gesicki et al. (2003). Together with the 33 new 
PNe, the full sample contains 101 objects (from 73 earlier ob- 
jects we removed three extragalactic PNe, and two PNe are re- 
analyzed with the new spectra). The full sample contains 23 
[WR], 21 wels and 57 non-emission-line central stars. 

Several objects were reclassified based on the work of 
Acker & Neiner (2003) and Gorny et al. (2004). The distinc- 
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Fig. 5. Stellar photo-ionization temperatures versus [WO]- 
[WC] subclass. 

tion between emission-line stars and non-emission-line stars 
depends on the depth of the available spectra. It is therefore 
possible that some of the non-emission-line stars would be clas- 
sified as wels stars with deeper spectra. The effect of detectabil- 
ity of the stellar lines on the classification is discussed by Gorny 
et al. (2004). 

4.1. Temperatures 

The subclass of the [WR] stars provides an indication of stel- 
lar temperature, but it can not be used to predict unique tem- 
peratures. This can for instance be seen for the earlier sub- 
classes (hotter stars) from the data in Acker & Neiner (2003). 
De Araujo et al. (2002) show that deep spectra are needed for 
accurate association of a [WR] subclass. In this paper we use 
the photoionization (equivalent black-body) temperature 7Vb, 
rather than the [WR] subclass. This parameter has the further 
advantage that it is also available for the non-[WR] stars. Fig.|5] 
shows the relation between photoionization temperature and 
[WR] ([WO] and [WC]) subclass. There is a good relation 
overall, but the change in temperature is smaller for the later 
subclasses (8-1 1). One discrepant object is M2-43 (27.6+04.2) 
where the temperature is much higher than the subclass ([WC 
7-8]) would suggest. The photo-ionization model fits the line 
ratios reasonably well, and the strong [Oin] line with weak 
[N n] require a high temperature (Acker et al. 2002). Such high 
temperature was also obtained by non-LTE model analysis of 
Leuenhagen & Hamann (1998). This object shows a dual dust 
composition (Gorny et al. 2001). The [WC] subclass should be 
checked. 

4.2. Expansion velocities 

For the whole sample of 101 PNe, the average (and the median) 
value of the mass-averaged expansion velocities is 22kms~ 1 
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with standard error on the mean of 0.67km s . The average V av 
for [ WR] PNe is 25 + 2 km s" 1 while for wels it is 22 ± 1 km s" 1 . 
The 1-cr difference in expansion velocity is much smaller than 
found in Pena et al. (2003a), who find an average expansion 
velocity for non-WR PNe of 21 kms~' (virtually the same as 
ours) but for [WR] PNe 36kms~ 1 . The difference can be ex- 
plained by their method of estimating the velocity: they use 
the half-width at half maximum of the nebular lines and this 
is affected by the turbulence in the velocity fields (Acker et al. 
2002). 

It is natural to explain the higher than average expansion ve- 
locity of [WR] PNe in terms of stronger stellar wind acting on 
the swept-up shell. Mellema & Lundqvist (2002) studied such 
interactions and confirmed the above tendency providing an ex- 
ample of [WR] PN expansion of 21 km s -1 versus 17 km s _1 for 
non-[WR] PN. This difference is smaller than obtained by Pena 
et al. (2003a), the values are rather comparable to ours. 

Pena et al. (2003b) argue that the expansion velocities for 
evolved [WR] PNe are larger than for younger objects. In that 
paper the expansion velocity is defined as half- width at 1/10 of 
the maximum intensity, a measure introduced by Dopita et al. 
(1988). At this intensity level the emission lines are strongly 
broadened by a turbulent component common for [WR] PNe 
so it is no surprise that our values (corrected for the turbulence) 
are smaller. Pena et al. (2003b) arrive at their conclusion by 
plotting stellar temperatures versus expansion velocities. In our 
sample, no such effect is visible (Fig.[6] left panel), except that 
the three highest expansion velocities in the sample are from 
[WR] stars. The V exp and 7Vb values are spread over the whole 
range. However, there is a trend for the [WR] stars to dominate 
towards the lower-right corner in the figure, i.e. lower values of 
Tb-b/Vexp- This trend looks very similar to that from Fig. 3 of 
Pena et al. (2003b). However their [WR] PNe clustering in the 
upper-right corner is not so similar. 

Schonberner et al. (2005a) discuss their observations of 13 
PNe and concluded that expansion speed increases with the 
r e fr of the central star, i.e. with the nebular age. Such a situ- 



ation cannot be explained with simple models which assume 
an initial nebular density distribution described by a power law 
p cc r~ a with a fixed index a. Their conclusion is that the 
power-law index increases systematically with T e g (see their 
Fig. 12). We plot their model sequences for four values of the 
power law index in the left panel of Fig. [6] The indexes are 
from left to right: 2.00, 2.5, 3.0 and 3.25. Our data for 101 PNe 
overlap with the lines (the velocities are not defined exactly 
in the same way). We do not confirm the Schonberner et al. 
(2005a) conclusion, but the range of power law indices which 
they derive provides a good fit to the spread seen in our sam- 
ple. The trend in Fig.[6]that [WR] stars dominate the distribu- 
tion towards lower values of 7Vb/ V exp , in combination with the 
slope of the model tracks, suggests that some [WR] PNe may 
have steeper initial density distributions than non-emission-line 
PNe. 

4.3. Dynamical ages 

The expansion velocity and radius of the nebula can be used 
to derive the age of the nebula. One can use the outermost ra- 
dius and the maximum expansion velocity, but this is affected 
by the acceleration caused by the overpressure in the ionized 
region. Instead we apply the mass-averaged velocity, and use 
for the radius 0.8 of the outer radius, roughly corresponding 
to the mass-averaged radius. To account for the acceleration of 
the nebula we average between the current velocity and that of 
the original outflow velocities on the Asymptotic Giant Branch. 
The procedure is described in Gesicki et al. (2003). We made a 
simplifying assumption about the AGB outflow velocity setting 
it equal to lOkms -1 for all objects. 

The dynamical ages can directly be compared between the 
different groups of stars. However, they may not be equal to the 
true ages of the nebulae, due to the various assumptions made. 
This problem was discussed together with non-LTE analysis 
of some central stars of PNe. Rauch et al. (1994) obtained for 
the object K 1-27 that the dynamical age was more than one 
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order of magnitude smaller than the evolutionary age. In an- 
other paper Rauch et al. (1999) analyzed four more PNe and 
obtained dynamical ages a few times larger than evolutionary 
ages. The nebular and stellar ages were compared on a much 
bigger sample by McCarthy et al. (1990). They obtained that 
generally the empirical ages were larger than the evolutionary 
ones however their method for estimating dynamical ages was 
recently strongly criticized by Schonberner et al. (2005b). The 
contemporary situation is far from being conclusive. 

In the recent article we found a nice support for our method 
for estimating nebular ages. Schonberner et al. (2005b) show 
that both radii and velocities increase gradually with time what 
justifies our averaging between original and current velocity. 
For ionization bounded models (dominant in our sample) the 
ages derived from ions [O in] and [N n] diverge in opposite di- 
rections from true ages what supports our multi-line analysis 
where such effects compensate. For density bounded PNe the 
situation is even better since both lines result in the same age 
equal to the evolutionary one. 

The dynamical age in combination with the stellar temper- 
ature gives the rate at which the central star is increasing in 
temperature. This rate is very sensitive to the core mass of the 
star, with more massive stars evolving much faster. To assign 
a core mass, we interpolate between published tracks (see e.g 
Gorny et al. 1997, Frankowski 2003). For a given dynamical 
age and stellar temperature we obtain the core mass (and lu- 
minosity). The dynamical ages and core masses are listed in 
Tabled As for the ages, the mass determinations are internally 
consistent between different objects, but may show systematic 
offsets with respect to the true masses. As a check of the con- 
sistency of the procedure the luminosities can be compared. In 
Gesicki et al. (2003) we obtained that the luminosities from 
photoionization modelling are in most cases smaller than those 
interpolated from evolutionary tracks. However this can be in- 
terpreted as leaking of the nebulae: the PN does not intercept 
and transform all stellar ionizing flux. This points to asymme- 
try or clumping. 

4.4. Turbulence 

For a number of nebulae, satisfactory fits to the velocity fields 
could not be obtained. In these cases we used enhanced tur- 
bulence, above what is expected from thermal broadening and 
instrumental resolution. Indicative of turbulence is that lines 
of all ions show Gaussian shapes of comparable width: ther- 
mal broadening gives much wider lines for the lightest element 
(hydrogen). Comparing the He with metal lines is here cru- 
cial because differences in thermal broadening between differ- 
ent metal ions are small. Therefore we searched for turbulent 
solutions only for those PNe with Ha and at least one metal line 
available. Turbulence was indicated for 30 of the 101 nebulae 
in the full sample. 

Gesicki & Acker (1996) and Acker et al. (2002) find that 
[WR]-type PNe are characterized by strong turbulent motions 
while non-emission-line PNe are not. This general tendency is 
confirmed in the present analysis, but there are some excep- 
tions. 



Among the wels we find both turbulent and non-turbulent 
objects. Our full sample contains 21 wels of which 5 show tur- 
bulence and 16 do not. Of these 21 objects, 10 come from 
the new observations described here (3 turbulent and 7 non- 
turbulent). 

We also find well-determined, turbulent solutions for 4 PNe 
not classified as [WR] nor wels (M2-16, M3-23, He2-128, 
He 2-129). Not all [WR]-type PNe in the new sample show tur- 
bulent motions: M 1-27 and M 1-37 have good solutions, with- 
out requiring turbulence. The two non-turbulent [WR] PNe are 
both classified by Gorny et al. (2004) as '[WC11]?'. They are 
the coolest objects in our sample. However, the strong corre- 
lation between emission-line stars and turbulence remains: the 
fraction showing turbulence for non-emission-line stars, wels 
and [WR] stars is 7%, 24% and 91%, respectively. 

Uncertainties in the [WR] identifications may affect these 
fractions. The turbulent non-WR PNe could have emission-line 
central stars. This is unlikely for M 2-16, a well-studied, highly 
metal-rich PN (Cuisinier et al. 2000). He 2-128 and He 2-129 
were studied spectroscopically by Dopita & Hua (1997): the 
non-WR identification cannot be attributed to a lack of obser- 
vations, but faint stellar lines may have been missed. Some ob- 
jects may have too few nebular lines to be able to detect turbu- 
lence: a velocity gradient broadens a single line profile similar 
to turbulence - we only accepted turbulent solutions where they 
fitted clearly better the line shapes than non-turbulent ones. 
Deeper spectra would help here. 

The onset of the turbulence and the [WR] phenomenon may 
not be completely simultaneous. However, we find no evidence 
that the turbulence increases for more evolved objects: Fig.|6] 
(right panel) shows no relation between turbulence and stellar 
temperature. 

In trying to understand why the presence of turbulence in 
wels PNe is so mixed, let us consider a possibility that a PN can 
switch between turbulent and non-turbulent regimes during its 
lifetime. The transition time for this should be comparable to 
the shock travel time across the nebula. Assuming a radius of 
0.1 pc and isothermal sound speed of 10kms _I we obtain a 
time of 10 4 years. Perinotto et al. (2004) stated that the shock 
front velocity can be 3—4 times faster than the isothermal sound 
speed: this reduces the time to a few thousands years. This is 
comparable to the ages of our PNe. In our sample we have a 
couple of (compact) objects with an age of about 1000 years 
which have already developed turbulence. Acker et al. (2002) 
proposed that turbulence in PNe is triggered or enhanced by 
[WR] stellar wind inhomogeneities. The nebulae surrounding 
very active [WR] cores are turbulent. The less active wels with 
much weaker winds, if they could change their wind charac- 
teristics, they would develop turbulence, or exhibit decaying 
turbulence, over the lifetime of the PN. Thus, the mixed oc- 
currence of turbulence in wels is compatible with a transitory 
status, not unlikely a [WR]-star progenitor, provided the life 
time of the phase does not exceed ~ 10 3 yr. 
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4.5. IRAS colours 

Planetary nebulae are strong infrared emitters, due to their 
heated dust. The dust colours follow the evolution: as the neb- 
ula expands, the dust cools and the dust colours redden. Zijlstra 
(2001) has shown that [WR] stars are stronger infrared emitters 
than other PNe, with a tendency for bluer colours. The bright- 
est infrared PNe are the IR-[WR] stars, a subgroup of the cooler 
[WC] stars. 

Fig.0shows the colour-colour diagram for the present sam- 
ple. Different symbols distinguish the [WR] stars (filled cir- 
cles), the wels (open circles) and the non-emission-line stars 
(pluses). In addition, overplotted star symbols indicate turbu- 
lent nebulae. We excluded the confused source M2-6 where 
the IRAS colours appear to be those of a nearby AGB star. 
The boxes are source classification regions of Zijlstra et al. 
(2001). Planetary nebulae fall in box II. The cooler [WC] stars 
are found mainly towards the left in the diagram, in three cases 
with colours like those of (young) post-AGB stars. (The three 
objects are BD +30 3639, M 4- 18 and He 2-142.) The wels are 
shifted to redder colours, and their distribution is similar to 
those of normal PNe. Note that the cooler non-emission-line 
stars in our sample have 60/25-micron colours similar to the 



cool [WC] PNe, but their 25/12-micron colours are very differ- 
ent. 

We constructed dust models following the nebular expan- 
sion occurring as the star heats up. We used a r~ 2 wind, with 
the model and recipe of Siebenmorgen et al. (1994). The model 
calculates the IRAS fluxes at any particular time. The results 
are shown in Fig.[7] for a carbon rich (long dashes) and a oxy- 
gen rich (short dashes) nebular dust. The evolutionary progres- 
sion from cool stars with bluer IRAS colours to hot stars with 
(more evolved) redder IRAS colours is visible for the [WR] 
stars , but not for the non-emission-line stars and the (less pop- 
ulated) wels group (Zijlstra 2001). 

Detection statistics show the differences between the 
groups. Out of the 101 stars, 60 have three good IRAS detec- 
tions. The detection rates are 28 out of 57 non-emission-line 
stars, 12 out of 21 wels and 20 out of 23 [WR] stars: respec- 
tively 49%, 57% and 87%. This shows that the [WR] stars are 
on average much stronger IRAS emitters. If we only consider 
25 and 60 micron, the detection rates become 89%, 86% and 
100%, i.e. the difference is rather reduced. The numbers indi- 
cate that [WR] stars are much brighter at 12-micron, but less so 
at the other bands. 

The wels have the same detection rate as other non-[WR] 
stars. However, all 5 wels with turbulence have three-band de- 
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distribution for hydrogen burning 0.605M o track of Blocker 
(1995). The predicted number of objects (per logarithmic tem- 
perature bin) gradually increases towards higher temperatures, 
as expected from the increasing width of the bins and from 
slowing the evolution. 

When compared to the observational data we see that no 
single group of PNe produces a histogram similar to the mod- 
els. The solid line in the bottom panel of Fig.|8]shows that the 
sum of all classes produce a histogram most similar to the the- 
oretically expected distribution. Thus, the full sample is repre- 
sentative for a uniformly sampled Blocker track. The flat dis- 
tribution of the emission-line stars implies that this group is 
biassed towards lower temperatures, relative to standard evolu- 
tion. The wels as a separate group show a strong temperature 
preference and cannot represent a full evolutionary track. 

It is possible that our observations misrepresent the oldest 
PNe with the hottest cores. Our sample is mainly restricted to 
the spherical PNe with an ionization boundary while for old 
PNe the ionization front breaks through and may cause frag- 
mentation of a nebula making it not suitable for our analysis. 
We also cannot exclude a systematic error of our 7Vb values 
which may shift the maximum of the observed distribution, and 
this error may also depend on temperature. 



5.2. Dynamical ages of the PNe 



tections; their colours are similar to the [WO] stars. All wels 
without 12-micron detection are non-turbulent objects. 

5. Evolutionary parameters 

5.1. Temperature distributions 

The temperature distribution of the three groups of stars, de- 
fined as in section l4~Tl is shown in Fig.|H] The [WR] stars show 
some clustering at low and high temperatures with a gap in be- 
tween (see upper panel). This gap has been noted before long 
time (see e.g. Gorny 2001). The wels stars, in contrast, peak 
precisely in this gap. In contrast, the non-emission-line stars 
show a peak in their distribution at high temperature (see cen- 
tral panel). 

Summing the two groups of emission-line stars (cen- 
tral panel), we find a flat distribution (in logarithmic bins!). 
Summing all three groups together (bottom panel, drawn line) 
shows an increasing number of objects with log T. 

We compare the temperature distribution of the various 
PNe samples with that predicted by the evolutionary tracks of 
Bloocker (1995). We interpolate the tracks to constant temper- 
ature steps, in order to calculate the predicted number of stars 
per temperature range (proportional to the time spent in the 
range). The model tracks have to first order a fairly constant 
temperature increase, dT/dt. In logarithmic bins, this means 
that the number of stars per bin increases with temperature. To 
be consistent with our model analysis we only use the the hori- 
zontal part of the evolutionary track, before the knee in the H-R 
diagram. Fig.[8](bottom panel, dotted line) shows the predicted 



Evolutionary tracks should connect objects along sequences of 
increasing dynamical age and increasing stellar temperature, 
until the stars reach the knee in the HR diagram, after which the 
stellar temperature slowly decreases. Helium-burning tracks 
generally also follows this sequence, but the stars evolve a few 
times slower (factor of 3 according to Blocker 1995); the details 
depend on the precise timing of the thermal pulse preceding the 
switch to helium burning. 

Fig. [9] plots the dynamical age versus the temperature for 
the stars in our sample. The interpolated (hydrogen-burning) 
evolutionary tracks are overplotted with a solid line. As in 
Fig.0different symbols distinguish types of stars and indicate 
turbulence. 

The figure shows that some regions in this age-temperature 
plane are dominated by certain types of objects. The [WR] stars 
are mainly found in the narrow region between the 0.605 and 
0.625 M tracks, with a few further objects at low temperatures 
and a range of ages. The wels stars are spread mostly out at 
intermediate temperatures. The non-emission-line stars tend to 
avoid the low temperature region, but cover a region in the top- 
right corner which the other stars tend to avoid. 

The statistics of our full sample of 101 PNe shows that the 
average dynamical age of non-emission-line central stars and 
of wels is 3200 years while for [WR] -stars it is 2800 years, 
with standard error on the mean of about 300 years in both 
cases. This difference is not significant and the result may be 
contaminated by a selection effect because it is easier to detect 
a [WR] star at lower temperatures (see Gorny et al. 2004). 



K.Gesicki et al.: Planetary nebulae with emission-line central stars 



11 




age [1000 yrs] 



Fig. 9. The age-temperature diagram. Filled circles indicate [WR] stars, open circles are wels and pluses indicate non-emission- 
line stars. Turbulent nebulae are indicated by the star symbols. Solid lines show Schonberner's H-burning evolution. Labeled are 
two individual PNe discussed in Sect.6.2. 



5.3. Non-emission-line objects 

In Fig.|9] the non-emission line objects are mixed with the 
emission-line ones except of the upper right area. This is the 
place where the cooling parts of evolutionary tracks of mas- 
sive stars cross with heating parts of evolutionary tracks of low 
mass stars. The hot (10 5 K), non-emission-line stars with old 
nebulae could be high-mass stars located on the cooling track. 
The alternative interpretation is that they have lower-mass cen- 
tral stars. In either case, the absence of emission-line stars in 
this region is obvious. Either a minimum mass is required for 
a hot [WO] star of about 0.61 M , or the stellar wind ceases 
when the star reaches the cooling track. The latter is plausi- 
ble since the luminosity drops quickly by a factor of 100 at 
this time. This fact is well known for hot star winds (see e.g. 
Owocki 1994, Owocki & Gayley 1995). Koesterke & Werner 
(1999) show that PG1159 stars have much weaker winds than 
wels stars, and suggest that the wind mass-loss rate decreases 
rapidly when the star enters the cooling track. 

5.4. Nebular evolution 

Schonberner et al. (2005a) conclude that during the evolution 
the PN material always enters the 'champagne phase' where 
the ionization front breaks through and the PN becomes den- 



sity bounded. This disagrees with our result that most of the 
101 PNe are ionization bounded despite their age and core T e g . 
In our sample the presence of the ionization boundary was de- 
duced from the strengths of the low-ionization low-excitation 
lines. A possible explanation of this disagreement can be found 
in Gesicki et al. (2003) where we compared the model lumi- 
nosities with those obtained from evolutionary tracks and sug- 
gested that 'PNe are leaking'. This indicates that PNe are a 
mixture of regions transparent and opaque for ionizing radi- 
ation. The spherical Schonberner's and Torun models cannot 
investigate asymmetric or clumpy nebulae. However there is 
another explanation - the ionization boundary is not so ex- 
ceptional. The evolutionary sequence No. 4 of Perinotto et al. 
(2004), with high AGB mass loss rate and low AGB wind ve- 
locity, never becomes optically thin and can serve here as a hint. 
It seems to be promising to follow this idea and to search for 
other hydrodynamical models that remain ionization bounded 
for a significant part of PN evolution. In Gesicki & Zijlstra 
(2003) and Zijlstra et al. (2005 in preparation) we presented 
models supported by HST images showing very small central 
cavities of the nebulae. This finding supports the above idea of 
a low AGB wind velocity, at least at the end of AGB and for 
some PNe. Finding an appropriate hydrodynamic model might 
constrain the still unknown AGB mass loss history. 
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6. The hydrogen-poor sequences 

6.1. The [WR] stars 

A global nebular quantity, turbulence, is predominant among 
the [WR] stars. Mellema (2001) links the metals-enriched 
[WR] winds with long lasting momentum-driven phase which 
drives the turbulence. This confirms that they form a distinct 
group, and it can be interpreted as evidence for a separate dif- 
ferent evolutionary channel. Fig.0and Fig.|5]show that there 
is a marked uniformity among the majority of the [WR] stars. 
Fig.|5]suggests this is related to the core mass of the star, where 
stars of a relatively narrow range in core mass account for the 
majority of the [WR] stars. Caution should be taken when quot- 
ing precise masses: the location in the age-temperature diagram 
depends on modeling assumptions as discussed above. 

The [WR] stars are in general hydrogen-poor, and are con- 
sidered helium-burning central stars. These stars will show a 
slower temperature increase than do hydrogen-burning tracks 
of the same mass, although the detailed evolution is compli- 
cated. The distribution of [WR] stars in Fig.|9]is not inconsis- 
tent with such an inclined track, but strong conclusions are not 
warranted. 

There are three [WC]-PNe located near the 0.565 M track: 
from the left M 4- 1 8, He 2-99 and NGC 40. All are old nebulae 
with cool stars. The first is determined as old because of small 
expansion velocity, the other two PNe are large («0.15pc) as 
expected for old objects. Their three central stars were analyzed 
by Leuenhagen et al. (1996). All have long been suspected to 
be born-again PNe, where the central stars have suffered a Very 
Late Thermal Pulse (VLTP), leading to a rejuvenated star in- 
side an old nebula (e.g. Hajduk et al. 2005). Helium-burning 
stars are not expected to experience another VLTP. The stars 
along the VLTP sequence therefore represent a separate evolu- 
tionary channel from the majority of the [WR] stars. The dif- 
ference may largely be one of timing, whether the last thermal 
pulse takes place along the horizontal sequence (LTP), or on the 
cooling track (VLTP). The only non-emission-line star in this 
region is the halo object PN G 359.2-33.5, which could well be 
a genuine low-mass object (McCarthy et al. 1991). The wels in 
this region is He 2-108 (its classification is uncertain). 

There is a group of 7 [WC] stars which combine very com- 
pact nebulae (low ages) with very cool central stars. This group 
has temperatures in a narrow range around log 7Vb = 4.7, well 
separated from hotter [WR] objects (see Fig.[5]and Fig.|9j. The 
IR-[WC] stars and the mixed-chemistry objects are found here 
(Zijlstra 2001). It is not clear that these belong to either of the 
two groups discussed above. The mixed chemistry has been as- 
sociated with binary evolution (Zijlstra et al. 2001, de Marco 
et al. 2004), and in one case (CPD-56 8032 — Cohen et al. 
1999, de Marco et al. 2002), there are strong indications for bi- 
narity. The precise role of a binary companion is not clear, but 
it may involve retaining old ejecta in a disk, triggering accre- 
tion and possibly providing a source of ionization. They could 
have similar origins to the R Cor Bor stars (see the review of 
Clayton 1996), although the very different circumstellar envi- 
ronments suggest the relation is not an evolutionary one. 



This group of cool, compact objects could in principle 
have descendants among the hot [WO] stars, which have more 
evolved nebulae. However, the mixed chemistry (i.e. evidences 
for both a neutral oxygen-rich and an inner carbon-rich re- 
gions) and unique IRAS colours argues against this. But there 
are no obvious descendants among the non-emission-line stars 
either, leaving the evolution of this group unexplained. If the 
stars are accreting, the extended, cool photosphere could hide 
a much hotter underlying star: this would explain the lack of 
intermediate-temperature related stars, but still leaves the ques- 
tion of descendants open. 

6.2. The gap in [WR] temperature distribution 

The lack of PNe cores of [WC 5-7] type is well known 
(Crowther et al. 1998, Gorny 2001). This feature can also 
be seen in our temperature distributions presented in Fig.|5] 
and Fig.|5] One of the two [WR] star filling this gap is PB 8 
(PNG 292.4+04.1) of [WC5-6] type (Acker & Neiner 2003). 
However by Parthasarathy et al. (1998) it was classified as 
wels and Mendez (1991) found it as H-rich. In Fig.Q this 
object is positioned among wels stars. Therefore we propose 
that PB 8 is a wels and no longer fills the [WR] temperature 
gap. Another object, the only known [WC 6] type - M 1-25 
(PNG 004.9+04.9) fits the centre of the gap. The Torun mod- 
els resulted in a stellar temperature log 7Vb = 4.6 (Gesicki & 
Acker 1996). However much higher temperatures have been 
published for the Ml-25 central star, reaching log 7Vb = 
4.9 (Samland et al. 1992) based on the helium ionization. In 
Fig. this object is placed amongst other regular [WR] ob- 
jects. Therefore we propose that M 1-25 is a [WC 6] object (we 
don't question the WR classification) but its temperature may 
be above the discussed gap. 

The uncertain status of the interlopers strengthens the case 
that the gap in the temperature distribution of [WR] stars is 
real and quite prominent: the region 4.55 < log 7Vb < 4.8 
appears devoid of [WR] stars. This supports the proposed sep- 
arate evolutionary sequences of cool and hot [WR] stars. The 
same conclusion was reached by Crowther et al. (1998) who 
applied a unified classification scheme to massive WR and low 
mass [WR] stars and realized that the mentioned above gap is 
not visible for WR objects. 

6.3. The wels 

The wels appear to be a mixed zoo. We tentatively identify three 
groups. First, a group which shows similarities in distribution 
in Fig.|5]to the [WR] stars. These include one object among 
the VLTP sequence, and approximately 5 objects which form 
a low-temperature extension to the hot [WO] stars and may 
be their progenitors. Second, around 6 wels are found along 
the high-mass track (> 0.625 M ). These we suggest are ob- 
jects where the high stellar luminosity drives a stellar wind- 
they may be H-rich and not related to the [WR] stars. Finally, a 
group along a low-mass track (M « 0.6 M Q ). These could rep- 
resent failed [WR] stars, where the luminosity is insufficient to 
drive a full WR-type wind. 
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The wels show on average a shift in location in Fig.|6] they 
fall in a region consistent with a lower power index than do 
[WR] stars. The nebulae around [WR] stars tend to be more 
centrally condensed than those of wels stars. This points to a 
different mass-loss history and is consistent with the sugges- 
tion that some fraction, possibly the majority, of wels are not 
evolutionary related to the [WR] stars. 

The small group which may be [WO] progenitors suggests 
a sequence wels— [WO], reflecting increasing temperature. This 
could reflect a bi-stability jump, where the strength of a stellar 
wind suddenly changes when a change in ionization equilib- 
rium makes different lines act in the absorption of radiation. A 
similar situation is known among early-type supergiants (e.g. 
Vink et al. 1999, Tinkler & Lamers 2002). But the latter au- 
thors find that among PNe cores, the wind suddenly decreases 
by a factor of up to 100, precisely at the temperature where 
we find the w/s-[WO] transition. They attribute this to a CNO 
bi-stability jump. Their analysis used stars classified as O(H) 
only, and could be biassed if hotter stars with high M appear as 
[WO] stars. 

For the other stars, the driving factor differentiating wels 
and [WR] stars appears to be luminosity, 

6.4. The PG11 59 stars 

The hot PG1 159 stars are generally considered the most natural 
descendants of [WR] stars (Werner 2001). They are not a uni- 
form group: Werner et al. (1996) refer to them as the 'PG1 159 
zoo'. There is also a group of so-called '[WC]-PG1159 tran- 
sition objects' which are found between [WCL] and [WCE] 
groups in the log g - log T e g plane (Werner 2001). Because our 
results suggest separate sequences for cooler [WC] and hotter 
[WO], one could envisage that cool [WC] stars evolve towards 
white dwarfs via a [WC]-PG1 159 phase, whilst hot [WO] stars 
do so via a PG1 159 phase. However, this remains speculative. 

It has been reported that the central star of the Longmore 4 
nebula changed from a PG1 159 type to a [WC 2-3], for a time 
of a few months (Werner et al. 1992). Its absorption spectrum 
transformed into an emission spectrum, with changes proceed- 
ing on a timescale of days. It was interpreted as strongly en- 
hanced mass loss, and may have been triggered by an accretion 
event. As yet it is the sole documented example of such impres- 
sive variability among PN cores, but even so it shows that the 
picture is a complicated one, where the chaotic picture of wels 
could be interpreted in terms of stellar wind activity. 

7. Conclusions 

We have deduced for 33 PNe the velocity field. Since this was 
done in most cases by considering two nebular lines we regard 
as robust the mass averaged value only. For six PNe the anal- 
ysis was based on three or more lines and is therefore more 
reliable. The discussion of details of velocity field like e.g. tur- 
bulence should be therefore treated with caution while the anal- 
ysis based on the mass averaged velocity (ages, masses) should 
be reliable. We combined this sample with earlier published 
data and discussed the full sample of 101 PNe. 



We found no evidence for an increase in nebular expansion 
velocity with time. This correlation was checked against 7Vb 
which should correspond to stellar age and on the nebular dy- 
namical age i.e. on two independent characteristics. The [WR] 
stars show a tendency to low ratios of 7Vb/Vexp- Comparison 
with the Schonberner et al. (2005a) models suggests that these 
objects have more centrally condensed nebulae, indicative of 
increasing mass-loss rate with time on the AGB. Emission-line 
stars are not found on the cooling track. 

We find a correlation between IRAS 12-micron excess and 
the [WR] stars. This excess is not shown by the wels. The cause 
of 12-micron excess is typically a population of small grains. 
These may form in the wind, or form in the collision between 
the stellar wind and the nebula. Without exception, the wels do 
not show the 12-micron excess. 

We confirm the stong correlation between turbulence and 
[WR] stars. However, the correlation is not one to one: there are 
a few non-emission-line stars with turbulence, and a few [WR] 
stars without. The wels appear intermediate in occurrence of 
turbulence. 

The distribution in the age-temperature plane suggests that 
there are several groups of [WR] stars and wels. One group 
shows evolved nebulae with cool stars and likely originates 
from a Very Late Thermal Pulse. The hot [WO] stars show 
a narrow distribution corresponding to a small range in core 
mass. Some wels fall on the same sequence but at lower 
temperature-they may be the [WO] progenitors. Other wels 
are distributed on higher and on lower mass tracks, which lack 
[WR] stars. Finally, a group of cool [WC] stars with IR-bright 
compact nebulae have no counterpart among any other group 
of PNe. 

The relation between the [WR] stars and the wels is there- 
fore a complicated one, with neither group being uniform. For 
the majority of the wels, there is evidence against them being 
evolutionarily related to the [WR] stars. There is a noticeable 
apparent difference in core mass and dust properties and neb- 
ular turbulence. A small group (5 objects) of wels have char- 
acteristics consistent with [WO] star progenitors. Apart from 
this wels-[WO] sequence, the wels contain stars at both the 
high and low luminosity end of the distribution. The former 
are likely stars where the strong radiation field drives the wind, 
but are not necessarily H-poor. The latter may be H-poor, failed 
[WR] stars where the luminosity is too low to maintain a strong 
[WR] wind. For the [WO] progenitors, the determining factor 
appears to be stellar temperature: the particular ionization bal- 
ance may reduce the efficiency of the radiation-driven wind. 

The determined discontinuity in temperatures of [WR]-type 
stars is in favor of the suggested separate evolutionary channels 
of cool and hot [WR] objects. 
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